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ABSTRACT 

We demonstrate for the first time that gaseous halos of disk galaxies can 
play a vital role in recycling metal-rich gas ejected from the bulges and thus 
in promoting chemical evolution of disks. Our numerical simulations show that 
metal-rich stellar winds from bulges in disk galaxies can be accreted onto the 
thin disks owing to hydrodynamical interaction between the gaseous ejecta and 
the gaseous halos, if the mean densities of the halos (phg) are as high as 10"^ 
cm~^. The total amount of gas that is ejected from a bulge through a stellar wind 
and then accreted onto the disk depends mainly on phg and the initial velocity 
of the stellar wind. About ~ 1% of gaseous ejecta from bulges in disk galaxies 
of scale length can be accreted onto disks around R ~ 2.5ad for a reasonable 
set of model parameters. We discuss these results in the context of the origin 
of the surprisingly high metallicities of the solar neighborhood disk stars in the 
Galaxy. We also discuss some implications of the present results in terms of 
chemical evolution of disk galaxies with possibly different phg in different galaxy 
environments. 



Subject headings: Galaxy: halo - Galaxy: bulge - Galaxy: disk - ISM: jets and 
outflows - galaxies:evolution - 
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1. Introduction 

Gaseous infall onto galactic disks from their halos has long been considered to play vital 
roles in various aspects of galaxy formation and evolution, such as formation of SO galaxies 
from truncation of halo gas infall (e.g., Larson et al. 1980), maintenance of spiral arms in 
disk galaxies (e.g., Sellwood & Carlberg 1984), and as a solution to the so-called G-dwarf 
problem (e.g., Pagel 1997). The total masses and abundance properties of halo gas infalling 
onto disks and the timescales of gas infall have been considered to be key determinants 
which can control galaxy evolution processes (e.g., ScIIwood & Carlberg 1984), in particular, 
chemical evolution of disk galaxies (e.g., Ghiappini et al. 1997). Although it remains largely 
unclear both observationally and theoretically how halo gas can be accreted gradually onto 
disks, the gradual infall of halo gas is one of key ingredients in most chemical evolution 
models that successfully explain their chemical properties (e.g., Matteucci & Francois 1989). 

It is well established that there are significantly metal-rich disk stars with [Fe/H]~ +0.4 
in the solar neighborhood (e.g., Feltzing & Gustafsson 1998; Bensby et al. 2005). Recent 
observational and theoretical studies have suggested that the canonical models of the Galactic 
chemical evolution based on infall of metal-poor halo gas have a serious problem in explaining 
the observed presence of metal-rich stars with [Fe/H] = -1-0.2—0.4 (e.g., Tsujimoto 2007). 
One of possible ideas to solve this problem is that the Galactic disk experienced accretion of 
metal-rich gas which originated from the bulge: the observed rather metal-rich disk stars in 
the solar neighborhood were formed from gas ejected from the bulge (Tsujimoto 2007). It 
is, however, totally unclear whether the gas ejected from the Galactic bulge can really reach 
the solar-neighborhood under hydrodynamical infiuences of the gaseous halo owing to lack 
of extensive numerical studies on dynamical fates of the ejected bulge gas. 

The purpose of this Letter is thus to show, for the first time, how gaseous halos of 
disk galaxies infiuence metal-rich stellar winds from the bulges through nuclear activity 
(i.e., starbursts and AGN). We numerically investigate whether and how stellar winds from 
bulges (or stellar ejecta due to supernova feedback) can be accreted onto the disks after 
hydrodynamical interaction between the ejecta and the gaseous halos. This investigation is 
important, not only because recent observational studies have found evidence of large-scale 
outflow from bulges in the Galaxy and other disk galaxies (e.g., Bland-Hawthorn & Cohen 
2003; Matthews & Grijs 2004; Veilleux et al. 2005), but also because the physical role of 
galactic gaseous halos in chemical evolution of disk galaxies has not been so far investigated 
by numerical simulations. We mainly discuss how the total amount and distribution of gas 
ejected from a bulge and then accreted onto the disk depend on physical properties of the 
gaseous halo and the bulge ejecta. It should be stressed here that the bulge wind would need 
to have taken place in the Galaxy several Gyr ago to explain the possible ages of the super 
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metal-rich solar neighborhood stars (Tsujimoto 2007). 

2. The model 

We numerically investigate the dynamical evolution of gas ejected from bulges in disk 
galaxies using our own GRAPE-SPH codes as we have done in previous studies (e.g., Bekki 
& Chiba 2006). A disk galaxy is represented by a dark matter halo, gaseous halo, stellar 
and gaseous disks, a stellar bulge, and gas initially within the bulge (referred to as "bulge 
gas" for convenience from now on). The dark matter halo has an "NFW" profile (Navarro 
et al. 1996) with total mass of M^m , scale length Vg, and c parameter of 10. The gaseous 
halo has mass Mhg and the same spatial distribution as the dark matter and is assumed to 
be initially in hydrostatic equilibrium. The initial gaseous temperature of a halo gas particle 
is therefore determined by the gas density, total mass, and gravitational potential at the 
location of the particle via Euler's equation for hydrostatic equilibrium (e.g., the equation 
lE-8 in Binney & Tremaine 1987). 

The bulge has a Hernquist profile (Hernquist 1990) with mass Mb and scale-length a\^. 
The disk has an exponential density profile with mass Mj, radius Rd, and scale length 
{— 0.2i?d). The stellar and gaseous disks have masses of Mg and Mg, respectively, and the 
same exponential density profile, and the gas mass fraction (Mg/M^) is a free parameter. 
In addition to the rotational velocity due to the gravitational field of disk, halo, bulge com- 
ponents, initial radial and azimuthal velocity dispersion are assigned to the disk component 
according to the epicyclic theory using Toomre's parameter (Binney & Tremaine 1987) Q 
= 1.5. We do not include the outflow of stellar winds from stars in the disk in the present 
models, because we think that the stellar winds with lower metallicities do not significantly 
infiuence chemical evolution of the disk in comparison with the metal-rich bulge wind. 

The initial distribution of the bulge gas with mass Mbg is assumed to be the same as 
that of the bulge and each bulge gas particle is assumed to be ejected from the bulge in 
the radial direction with the ejection velocity fcj. Also the particles are assumed to have 
initial rotational velocities of frot (around the z-a.xis) ranging from to a circular velocity 
at R — Sob. We show the results of the models with Vrot — 212 km s~^, because the results 
do not depend strongly on v^ot- All bulge gas particles are assumed to have the same initial 
Vej and I'rot in each simulation. The initial temperature of the bulge gas (Tbg) is assumed 
to be 10^ K, which corresponds to warm ionized gas observed in star burst galaxies such as 
M82 (e.g., Heckman et al. 1987). The models with very high temperatures of Tbg ~ 10^ K 
do not show significant gas accretion onto disks. 
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Since we mainly discuss the abundance properties of the Galactic disk, we chose the 
following values of the model parameters for disk galaxies: = lO^^M©, rg/od = 5, 

Md = 6 X lOi°M0, ad = 3.5 kpc, /g = 0.1, Mb = IO^^Mq, and Ob = 0.7 kpc. The most 
important parameters in the present study are the mean halo gas density p^g within SR^ 
determined from Mhg and Vcj for a given Mbg. Although we have investigated many models 
with different phg and v^j, we show the results of the standard model Ml with phg = 10~^ 
cm~^, = 500 km s~^, and Mbg = IO^Mq. Strong bulge winds can be more clearly seen in 
the models with Vej > 500 km s~^. Sembach et al. (2003) showed that the Galaxy gaseous 
halo is highly extended {R > 70 kpc) and low-density (phg < 10~^ — 10~^ cm~^), which means 
that the adopted phg is reasonable. The total number of particles used in a simulation ranges 
from 228805 to 258805 depending on /g and Mhg. 

Parameter values of the models (Ml to M7) discussed in this paper are shown in Table 
1. We do not discuss how rotation of gaseous halos influences the accretion processes of 
bulge gas onto disks, because we find that the influence is not significant. We show some 
results of comparative models in which fg — and phg = cm~^: these models enable us to 
grasp essential roles of gaseous halos in the accretion processes, though they are unrealistic. 
We will describe in detail the results of models not discussed in this paper (e.g., those with 
different Mbg and Vrot) in our future papers. We mainly investigate the final total mass 
of bulge gas particles settled around "the solar neighborhood" {R — Rq ~ 2.5ad) with 7 
kpc < i? < 10 kpc and \z\ < 1 kpc (Mace) for each model. We also investigate the mass 
ratio (/acc) of Mace to Mbg for each model. The final gaseous distributions at T = 0.45 Gyr 
correspond to those of the Galaxy about several Gyr ago in the present study. 



3. Results 

Fig. 1 shows how gas ejected from a bulge evolves with time during hydrodynamical 
interaction between the gas and the halo gas in a disk galaxy for the standard model Ml. 
Although a significant fraction of the gas once escapes from the bulge (T = 0.11 Gyr), most 
of the gas can be finally returned to the disk plane owing to gaseous pressure from the halo 
and the disk (T = 0.45 Gyr). The bulge gas particles initially in the central region of the 
bulge can interact so strongly with the disk gas particles immediately after the ejection that 
a thin disk can be formed from the bulge gas at i? < 3 kpc for a short time scale (T = 0.06 
Gyr). The bulge gas particles ejected from the bulge regions with larger \z\ can be later 
accreted onto the outer part of the disk and consequently can rotate around the center of 
the galaxy. The mean orbital eccentricity of bulge gas particles in the solar neighborhood is 
0.08 in this model, which means that the orbits become almost circular, because the particles 
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have acquired orbital angular momentum during hydrodynamical interaction with the disk 
gas. 

Fig. 2 shows that Mace increases with time owing to rapid accretion of the bulge gas and 
the mean accretion rate within 0.45 Gyr is 2.2 x 1O~^M0 yr""^ in the standard model. The 
mass-ratio (/acc) of Mace to Mbg is ~ 0.01, which means that a small fraction of the bulge gas 
can be accreted onto the solar neighborhood. Fig. 2 also shows that this result of small Mace 
is true for the comparative model M2 with reasonable pi^^, which suggests that the bulge gas 
can not so efficiently be accreted onto the disk outside the solar neighborhood irrespective 
of phg- The more rapid accretion and smaller Afacc in the model M2 is due largely to the 
stronger hydrodynamical interaction between the bulge and the halo gas. Fig. 3 shows that 
almost 80% of the bulge gas can be accreted onto the (inner) disk region with R <3 kpc in 
the standard model: this preferential accretion can be clearly seen also in the model M2. 

Fig. 4 compares the final distributions of bulge gas particles at T = 0.45 Gyr between 
different four models (M2 — 5). Only a compact thin disk with Mace = 1.4 x lO^Mo can be 
developed from the bulge gas in the model M3 without halo gas, which clearly demonstrates 
that hydrodynamical interaction between bulge and halo gas in a disk galaxy is crucial for 
the accretion of the bulge gas onto the outer part of the disk. Model M4 contain no disk and 
halo gas (i.e., /g = and p^g = cm^"^) and does not show a thin disk, which means that the 
presence of gas disks is also important for recycling of the bulge cjccta in disk galaxies. This 
result of M4 combined with that of Ml confirms that hydrodynamical interaction between 
bulge and disk gas is essential for forming a rotating gas disk composed of the bulge gas. 

The bulge gas can be pushed back to the inner disk more strongly and quickly by the 
halo gas owing to stronger hydrodynamical pressure of the gaseous halo against the bulge 
gas in the model M2. In this case the bulge gas can be accreted only onto the inner disk 
rather than the outer one in M2: only a very small fraction of the bulge gas can reach the 
solar neighborhood and outer gaseous streams and small gas clumps seen in the standard 
model are not seen in M2. The small gas clumps seen in the standard model might well be 
identified as high-velocity clouds, if they can form H I gas owing to efficient cooling. Less 
energetic gaseous flow from the bulge can be pushed back to the disk by the gaseous halo so 
that a thin disk with no small gas clumps is formed in the model M5. 

The dependences of /aec on model parameters are summarized as follows. Firstly, there 
is an optimum Vej which maximizes /aec for a reasonable phg = 10~^ cm~^. A larger fraction 
of bulge gas can escape from the disk for larger Vej so that a smaller fraction of the ejecta is 
finally accreted onto the solar neighborhood in the model M6 with Vgj = 1000 km s~^. On 

the other hand, a smaller fraction of the bulge gas can reach R = Rq owing to less energetic 
fiow so that a smaller fraction of the gas can be accreted onto the solar neighborhood in 
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model M5 with v^j — 200 km s ^. Thus, face is maximum for moderately energetic stellar 
winds such as Vej — 500 km s~^. 

Nevertheless, /acc can be higher for higher phg in the models with Vej — 1000 km s~^: 
in this case, higher densities of the halo gas are required for preventing the bulge gas from 
escaping the disk for more energetic stellar winds. This phg-dependence of face is totally 
different from that seen for Ve] = 500 km s~^, which implies that face depends on phg and Vgj 
in a complicated way. The model M4 with no halo/disk gas shows face (0.18) higher than 
that in model M3 with disk gas. Such high face in M4 is due to bulge gas particles that 
happen to be around the solar-neighborhood but are not rotating the disk: these particles 
soon return back to their original locations. It should be finally noted that the present results 
do not depend strongly on Vj-ot- 



4. Discussion and conclusions 

Our results indicate for the first time that gaseous halos of disk galaxies are very impor- 
tant for chemical evolution of the disks in the sense that they enable the galaxies to recycle 
the metal-rich stellar winds from the bulges and accelerate the chemical enrichment of galac- 
tic disks. The derived preferential accretion of the gaseous ejecta onto the inner disks will 
lead to a steepening of the abundance gradient after the action of strong central starbursts 
in the bulges followed by stellar winds. This steep abundance gradient is predicted to flatten 
with time owing to a subsequent chemical enrichment Tinder an accretion of low-metallicity 
infall from the halo (Tsujimoto et al. 2008), which is in accord with the observed flattening 
feature in radial metallicity gradients over the last several Gyr (Chen et al. 2003; Maciel et 
al. 2006, 2007). 

Our important flnding is that the metal-rich gaseous ejecta from the bulge can reach 
and chemically enrich the solar neighborhood through a one-time sporadic accretion event 
triggered by a starburst. This mechanism offers a promising channel to produce super-metal- 
rich stars with -1-0.2 <[Fe/H]< -1-0.4, the presence of which can not be accounted for by the 
conventional scheme of Galactic chemical evolution models (Tsujimoto 2007). An important 
caveat is that the observed very metal- rich stars with [Fe/H]~ -1-0.4 would need to be formed 
directly from the bulge gas without significant dilution with metal-poor halo and disk gas. 

Using the observed metallicity distribution function by Nordstrom et al. (2004), we 
estimate that the mass fraction of metal-rich stars with -1-0.2 <[Fe/H]< -1-0.4 is about 4% in 
the solar annulus. This means that the possible total stellar mass (Mmr) of the metal-rich 
stars for 8 kpc < -R < 9 kpc (i.e., the solar annulus) is about 3 x IO^Mq for the adopted 
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Galactic disk mass of 6 x lO^^M© and the exponential scale length of 3.5 kpc. The total 
mass of the metal-rich stars thus would give some constraints on (i) the total mass of the 
initial bulge wind (thus the strength of the central starburst) and (ii) to what degree the 
metal-rich bulge gas (e.g., [Fe/H]~+0.4) needs to be diluted via mixing of halo/disk gas to 
form less metal-rich stars with [Fe/H]~+0.2. 

Our study shows that about 1 % of gas ejected from bulges in disk galaxies can be 
accreted onto the disks at 7? 2.5ad for a reasonable set of model parameters in the present 
study. This suggests that chemical evolution of the outer parts {R ~ 3a^) of the disks can 
not be so strongly influenced by the accretion of the metal-rich gaseous cjccta. This also 
implies that if infalling metal-rich gas with [Fe/H]~-(-0.4 in the Galaxy can be mixed well 
with the gas with [Fe/H] ~ and a similar mass already present in the solar neighborhood 
and then form new stars from the mixed gas, the new stars would have [Fe/H]~-|-0.2: they 
however are unlikely to exhibit super-solar metallicities with [Fe/H]~ -1-0.4. However, if new 
stars can form preferentially from the metal-rich gas during the accretion of the gas onto the 
disk for some physical reasons, formation of new stars with super-solar metallicities would 
be possible. Our simulations imply that (i) dilution of the metal-rich bulge gas is likely to 
occur after the accretion of the gas onto the Galaxy and (ii) the dilution due to interaction 
with halo gas would be unlikely owing to rapid accretion of the bulge gas. 

Previous chemical evolution models of disk galaxies did not consider rapid accretion of 
metal-rich gas ejected from bulges onto disks (e.g., Laccy & Fall 1985; Prantzos & Aubcrt 
1995; Chiappini et al. 1997). The present study suggests that future more sophisticated 
chemical evolution models need to include possible influences of the rapid accretion on chem- 
ical evolution for disk galaxies, in particular, for those with bigger bulges which would have 
experienced stronger nuclear starbursts. The present study has not clarified (i) local (pc- 
scale) star formation processes during the accretion of metal-rich gas onto disks and (ii) 
metallicities of new stars formed from the star formation. We plan to investigate numeri- 
cally how mixing processes of metal-poor disk gas and metal-rich bulge ejecta within disks 
determine the metallicities of new stars formed during the accretion of the ejecta. 

We are grateful to the anonymous referee for valuable comments, which contribute to 
improve the present paper. K.B. acknowledges the financial support of the Australian Re- 
search Council throughout the course of this work. The numerical simulations reported here 
were carried out on GRAPE systems kindly made available by the Center for Computational 
Astrophysics (CfCA) at National Astronomical Observatory of Japan (NAOJ). 
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Table 1. Model parameters and brief results. 



model no 


Phg (cm 3) 


Wej (km s ^) 


/acc (xlO-2) 


comments 


Ml 


10-^ 


500 


1.03 


the standard model 


M2 


10-^ 


500 


0.39 




M3 





500 


0.02 


without halo gas 


M4 





500 


0.18 


without halo/disk gas (/g = 0) 


M5 


10-' 


200 


0.22 




M6 


10-5 


1000 


0.18 




M7 


10-^ 


1000 


0.55 
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Fig. 1. — Time evolution of the mass distribution in a disk galaxy projected onto the x-z 
plane for the standard model Ml with = 10~^ cm~'^. The time T in units of Gyr is 
given in the upper left corner for each panel. The mass distributions of disk stellar particles 
and halo gaseous ones are shown by blue and red, respectively, in the upper left panel for 
T = 0. For other panels, only the distribution of gas ejected from the bulge (i.e., "bulge gas 
particles" ) is shown by magenta so that the time evolution of the distribution can be more 
clearly seen. Note that most of the bulge gas particles can finally settle down to the thin 
disk after hydrodynamical interaction between the bulge and halo gas. 
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Fig. 2. — The time evolution of the total amount of bulge gas particles accreted onto the 
thin disk with 7 kpc < it! < 10 kpc (i.e., around the solar neighborhood, R — Rq) and with 
\z\ < 1 kpc in the standard model with phg = 10~^ cm~^ (solid) and the model M2 with 
Phg = 10"^ cm-3 (dotted). 
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Fig. 3. — The cumulative mass of bulge gas particles within R that are located within the 
thin disk (i.e., \z\ < 1 kpc) at T = 0.45 Gyr for the standard model (solid) and the model 
M2 (dotted). 
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Fig. 4. — Final distributions of bulge gas particles projected onto the x-z plane at T = 0.45 
Gyr for the model M3 without halo gas (upper left), M4 without halo/disk gas (upper right), 
M2 with phg = 10~^ cm-3 (lower left), and M5 with ^;ej = 200 km s'^ (lower right). 



